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(57) Abstract: The invention relates to a method of visualizing a three-dimensional volume image data set (1). A two-dimensional 
image (15) is then formed by projecting the volume image data on a projection plane in that for each pixel ( 1 6) there is determined an 
intensity value which corresponds to the maximum value or minimum value of the mean intensity along a projection path (17) which 
extends through the imaged volume and is associated with the relevant pixel. In order to determine the mean intensity, averaging 
is performed over a plurality of intensity values which neighbor one another along the projection path (17). In order to achieve 
effective noise suppression as well as adequate image contrast, averaging is performed over at least partly overlapping intervals (18 
to 22) along the projection path (17). The method in accordance with the invention offers the advantage that the user can choose 
the width of the averaging interval as an additional free parameter, so that image contrast and background noise can be influenced 
independently from one another to a significant extent. 
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The invention relates to a method of visualizing a three-dimensional volume 
image data set which contains a plurality of intensity values on a discrete spatial grid, in 
which method a two-dimensional image is formed by projecting the volume image data on a 
projection plane in that for each pixel there is detennined an intensity value which 
corresponds to the maximum value or minimum value of the mean intensity along a 
projection path which extends through the imaged volume and is associated with the relevant 
pixel, the mean intensity being calculated by averaging over a plurality of intensity values of 
the volume image data set which neighbor one another along the projection path. The 
invention also relates to an imaging diagnostic apparatus, notably a CT apparatus or an MR 
apparatus, and to a computer program for carrying out the method in accordance with the 
invention. 

Optimum visualization of image data is of crucial importance for medical 
applications, because such applications concern the formation of diagnostic images which are 
produced, for example, by means of computed tomography (CT) or by means of magnetic 
resonance tomography (MR). The object is to make the characteristic anatomical structures 
rnsxde a patient to be examined recognizable on the basis of a two-dimensional rendition. In 
medical imaging volume image data of a region of interest of the patient to be examined is 
reconstructed from the absorption of the X-rays (CT) or from the spin resonance signals 
(MR). The resultant image data set contains intensity values for each point of a discrete 
spatial grid. This data set, consisting of equidistant, so-called voxels, is further processed by 
means of suitable methods in the reconstruction unit of the relevant diagnostic imaging 
apparatus, that is, so as to visualize the anatomical structures, consisting of different types of 
tissue, on the basis of the different image intensities. 

The quality of the visualization is decisively dependent on the method used to 
transform the three-dimensional volume image data into a two-dimensional rendition. It is 
important that no crucial image information is lost by the reduction from three to two 
dimensions. 

There are various known methods of forming a two-dimensional image by 
projecting the volume image data on a projection plane. Each pixel of the projection image is 
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then associated with a rectilinear projection path which extends through the imaged volume 
and is oriented perpendicularly to the image plane. For each pixel an intensity value is 
determined by evaluation of the intensity variation along the relevant projection path. In the 
simplest case the intensity values of the volume image data set are averaged along each 
5 individual projection path. The intensity values are then summed and subsequently divided 
by the number of voxels situated along the relevant projection path. This approach offers the 
advantage that background noise which is inevitably present in real image data is effectively 
suppressed as a result of the averaging. This method, which is also referred to as the 
Collapsed View (CV) method, has the unacceptable side effect that the image contrast is 

10 reduced to such an extent that only objects exhibiting extreme differences in intensity can 
still be recognized in the resultant image. Depending on the relevant application, the image 
structures of interest are characterized by particularly high or particularly low intensity 
values. CT angiography, for example, utilizes the feet that the blood present in the vessels 
being examined yields a higher voxel intensity than the other organs of the imaged anatomy. 
.15 In these cases the use of the so-called Maximum (Minimum) Intensity Projection (MIP and 
mlP, respectively) method has become customary for generating two-dimensional projection 
images. The intensity value which corresponds to the maximum value or the mini?™ mi value 
of the voxel intensity along the projection path associated with the relevant pixel is then 
associated with each pixel of the two-dimensional projection image. The MIP or mBP method 

20 is very fast and reliably produces the desired image information upon reduction of the image 
data to two dimensions. However, it is also known that this method tends to give rise to 
artifacts which cannot be readily recognized in the resultant image. Furthermore, the MIP or 
mlP method emphasizes the background signal (noise) so that in given circumstances 
individual image objects will not be clearly visualized in the resultant image. 

25 In order to overcome said drawbacks, Thelissen proposes a so-called 

Maximum Average Projection (MAP) method in which averaging over a selectable number 
of voxels which succeed one another along the relevant projection path is performed prior to 
the determination of the maximum intensity value (Guillaume R.P. Thelissen, "Maximum 
Average Projection (MAP): A fast and robust way to visualize 3-dimensional MR data-sets", 

30 Proceedings of the ISMRM 1998, 2105). The known method at the same time effectively 
suppresses the background noise and the desired contrast enhancement is achieved as in the 
MIP method Since it is possible to select the number of successive voxels along the relevant 
projection path over which averaging takes place in accordance with the known method, the 
quality of the two-dimensional projection image is adjustable. For adaptation to the relevant 
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image noise the user can choose this value to be such that an as optimum as possible 
compromise is obtained between noise and image contrast 

The main drawback of the known method resides in the feet that as before a 
contradiction continues to exist between the desire for an as weak as possible image noise 
and an as high as possible image contrast. According to fee known method fee volume image 
data set is subdivided into mutually independent image slices prior to fee determination of fee 
maximum intensity values, fee thickness of said slices being dependent on fee number of 
voxels over which fee averaging is carried out The previously mentioned Collapsed View 
(CV) method is applied to a given extent to each of these image slices in order to calculate 
fee mean intensity values. These mean intensity values are compared wife one another along 
fee relevant projection path in order to determine fee maximum value for fee resultant 
projection image. The thicker fee image slices are, feat is, fee larger fee number of voxels 
over which averaging takes place, fee smaller fee number of mean intensity values evaluated 
during the projection will be. When a particularly strong noise suppression is desirable, 
averaging should take place over intervals which are as large as possible. However, mil has 
an adverse effect on fee image contrast because fee projection is then based on a 
correspondingly smaller number of mean intensity values. It is also particularly 
disadvantageous that in fee known method image information may be lost when fee 
structures imaged are situated at fee boundary between neighboring averaging intervals. This 
may give rise to image artefacts so that fee known method becomes rather problematic, that 
is, notably for diagnostic imaging. 

Considering fee foregoing, it is an object of fee present invention to improve 
fee known method to such an extent that an as effective as possible noise suppression takes 
place without giving rise, if possible, to loss of image information. 

This object is achieved on fee basis of a method of fee kind set forth in feat 
averaging is carried out over at least partly overlapping intervals along fee projection path. 

An advantage of fee method in accordance wife fee invention resides in fee 
feet feat because of their partial overlapping, fee averaging zones are larger than in fee 
known method, so feat overall fee noise suppression is more effective. A further advantage 
consists in feat as opposed to fee known method, fee number of mean intensity values 
evaluated during fee projection is not limited by fee selected width of fee averaging interval. 
Because fee averaging intervals overlap at least partly, fee number of intensity values over 
which fee averagmg operation is carried out can be selected independently of fee number of 
mean values. Because fee density of fee mean intensity values along fee projection path may 



WO 03/009227 PCT/EB02/03091 

4 

be significantly higher, the desired maximum or minimum intensity value will be determined 
with an accuracy which is significantly higher than that achieved in the known method. 
Moreover, no undesirable image artefacts occur, not even when the width of the averaging 
interval is larger than the image structures of interest 
5 In an advantageous further version of the method in accordance with the 

invention the averaging is carried out by convolution of the intensity variation along the 
projection path with a weighting function which has a selectable width which is larger than 
the distance in space between two neighboring data points of the volume image data set Such 
a convolution provides smoothing of the intensity variation along the projection path. The 

10 width of the smoothing function can then be selected at will as an independent parameter. 
The method in accordance with the invention can be carried out in such a manner that, 
despite the averaging, the overall number of intensity values on which the projection is based 
is maintained. When a suitable weighting function is chosen for the convolution, for example, 
a Gaussian function, the maxima or minima of the intensity variation will be maintained even 

1 5 when the averaging interval is significantly larger than the image structures to be imaged. 

It has been found that for a practical application of the method in accordance 
with the invention it is effective to cany out the averaging over intervals which overlap each 
time by one half. Using the same number of mean intensity values which are evaluated 
during the projection, the width of the averaging interval is then twice as large as in the 

20 known method. The noise suppression is thus significantly improved whereas the image 
contrast is not noticeably degraded. 

A CT apparatus or an MR apparatus or other imaging diagnostic apparatus as 
claimed in claim 4 is suitable for carrying out the method in accordance with the invention. 
The method in accordance with the invention can be advantageously implemented, that is, 

25 without necessitating special adaptations of the hardware, in conventional diagnostic 
apparatus in clinical use, since only the reconstruction unit is provided with a suitable 
program for the visualization in accordance with the invention. 

A computer program as disclosed in the claims 5 and 6 is suitable for this 
purpose. A computer program of this kind can be made available to the users of conventional 

30 imaging diagnostic apparatus on a suitable data carrier, for example, a disc or a CD-ROM, or 
can be offered for downloading via a public data network (the Internet). 

Embodiments of the invention will be described in detail hereinafter with 
reference to the Figures. Therein: 

Fig. 1 illustrates the projection method in accordance with the invention; 
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Fig. 2 shows the averaging with averaging intervals which overlap by one half 
and " * ' 

Fig. 3 shows a CT apparatus in accordance with the invention. 
Fig. 1 shows a three-dimensional volume image data set which is denoted by 
the reference numeral 1. The examination volume shown has an approximately cubic shape 
and contains a blood vessel 2 which is characterized by a particularly high image intensity 
The volume image data set 1 consists of twelve parallel slice images 3 to 14. Each of the slice 
images 3 to 14 consists of a two-dimensional matrix of intensity values. In order to visualize 
the volume image data set 1, a two-dimensional image 15 is formed by projecting the volume 
.mage data on Ihe plane of the projection image 15. For each pixel 1 6 there is determined an 
intensity value which corresponds to the maximum or minimum value of the mean intensity 
along a projection path 17 which extends through the imaged volume and is associated with 
the relevant pixel 1 6. The plane of the projection image 15 is oriented parallel to the slice 
images 3 to 14 in Fig. 1. The projection path 17 extends rectilinearly through the volume 
image data set and is oriented perpendicularly to the projection plane. For toe detennination 
of the intensity value for the pixel 16 first the mean intensities within the volume image data 
set 1 are calculated by averaging over a plurality of intensity values which neighbor one 
another along the projection path 17. Fig. 1 shows five averaging intervals 18 to 21 which 
partly overlap one another in conformity with the invention. For each of the averaging 
intervals 1 8 to 21 there is determined a mean intensity value. Thus, five mean values are 
compared with one another in order to find the maximum or minimum value of the mean 
intensity. For the interval 18 a mean intensity value is calculated for the corresponding point 
on the projection path 17 by averaging over the slice images 3, 4, 5 and 6. The mean intensity 
value for toe next mtervall^ ^ fa 

accordance with the invention the slices 5 and 6 are used for the determination of the mean ' 
value in the interval 18 and in the interval 19. The method in accordance with toe invention 
yields a planar image 22 of the blood vessel 2. 

Fig. 2 shows various possibilities for the averaging by means of averaging 
intervals which overlap each time by one half, that is, in the case of a volume image data set 
consxstmg often slices. In the upper diagram averaging takes place over the entire thickness 
23 of the examination zone. This approach corresponds to the previously described Collapsed 
View (CV) method. In the diagram shown therebelow averaging takes place over two 
intervals 24 and 25 which overlap in the central third part of the overall examination volume 
Three, four and nine mean intensity values, respectively, are determined in the diagrams 
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which are shown therebelow. In the case of nine averaging operations the width of the 
averaging interval 26 corresponds exactly to two image slices. In the diagram shown at the 
very bottom no averaging takes place. In that case the projection takes place in conformity 
with the customary MDP or mlP method. 

The CT apparatus as shown in Fig. 3 consists of a portal 27 on which an X-ray 
source 28 rotates about a patient table 29 and a patient 30 arranged thereon. On the portal 
there is also mounted a radiation detector 3 1 which faces the X-ray source 28 and also rotates 
about the patient 30. The radiation detector 3 1 detects the attenuation of an X-ray beam 
emitted by the X-ray source 28. The X-ray intensity is further processed and digitized by a 
measuring device 32. This data is then applied to a microcomputer 33 which reconstructs 
volume image data from the X-ray absorption data. The computer 33 is programmed in such 
a manner that in conformity with the invention it forms a two-dimensional image 1 5 from the 
reconstructed volume image data, said image being displayed on the display unit of the 
computer 33. 
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1 • A method of visualizing a three-dimensional volume image data set (1) which 

contains a plurality of intensity values on a discrete spatial grid, in which method a two- 
dimensional image (15) is formed by projecting the volume image data on a projection plane 
m that for each pixel (16) there is determined an intensity value which corresponds to the 
maximum value or minimum value of the mean intensity along a projection path (17) which 
* associated with the relevant pixel and extends through the imaged volume, the mean 
mtensity being calculated by averaging over a plurality of intensity values of the volume 
unage data set (1) which neighbor one another along the projection path, characterized in that 
the averagmg is carried out over at least partly overlapping intervals (18 to 22) along the 
projection path (17). 

2. A method as claimed in claim 1, characterized m that the averaging is carried 
out by convolution of the intensity variation along the projection path with a weighting 
function which has a selectable width which is larger than the distance in space between two 
neighboring data points of the volume image data set. 

3. A method as claimed in claim 1, characterized in that the intervals over which 
the averaging is performed overlap each time by one half. 

4. An imaging diagnostic apparatus, notably a CT apparatus or an MR apparatus 
for carrying out the method claimed in claim 1, which apparatus includes an imaging unit 
(27, 28, 31, 32) for the acquisition of coarse data of an object (30) to be examined and also 
mcludes a program-controlled reconstruction unit (33) which reconstructs volume image data 
from the coarse data, said volume image data consisting of a plurality of intensity values 
discrete spatial grid, and also forms a two-dimensional image (15) therefrom by projecting 
foe volume image data on a projection plane, characterized in that projection of foe volume 
image data is carried out by foe reconstruction unit (33) in such a manner that for each pixel 
(16) there is determined an intensity value which corresponds to foe maximum value or 
rninimum value of foe mean intensity along a projection path (17) which extends through foe 
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imaged volume and is associated with the relevant pixel, the mean intensity being calculated 
by averaging over at least partly overlapping intervals (18 to 22) along the projection path 
(17). 



5 5. A computer program for carrying out the method claimed in claim 1 , 

characterized in that it forms a two-dimensional image (15) from a three-dimensional volume 
image data set (1), consisting of a plurality of intensity values on a discrete spatial grid, by 
projecting the volume image data on a projection plane in that for each pixel (16) it 
determines an intensity value which corresponds to the maximum value or minimum value of 
10 the mean intensity along a projection path (1 7) which extends through the imaged volume 
and is associated with the relevant pixel, the mean intensity being calculated by averaging 
over at least partly overlapping intervals (1 8 to 22) along the projection path (17). 

6. A computer program as claimed in claim 5, characterized in that it performs 

15 the averaging by convolution of the intensity variation along the projection path with a 

weighting function which has a width which can be selected by the user and is larger than the 
distance in space between two neighboring data points of the volume image data set 
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